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1. INTRODUCTION 

The search for the scalar particles referred to as "Higgs bosons" [|l| is one of the major 
goals of present and future colliders. Such particles break the electroweak symmetry and 
are responsible for the masses of the fermions and bosons. The Standard Model (SM) 
predicts one neutral scalar (0°) while extensions of the SM often predict several scalars, 
both neutral and charged |^. Since 1989 the e^e~ collider LEP has searched for (fp. In 
the final run at energies around ^/s = 208 GeV two of the four experiments found excesses 
in the search for a Higgs boson in association with a Z boson |Q. Although it was agreed 
that a further run with about 200 pb~^ per experiment at a centre-of-mass energy of 
208.2 GeV would enable the four experiments to establish a 5a discovery (assuming the 
signal is genuine) f^ , the extended run was not approved and LEP was consequently shut 
down. The search for Higgs bosons will continue with Run II at the Tevatron which 
has a chance of confirming the existence of the Higgs boson in the mass range hinted at 
by LEP (?^ 115 GeV). This region is fairly problematic for the LHC and would require 
several years searching in the channel /i — ;> 77 to confirm such a light Higgs. A higher- 
energy e'^e" collider (NLC) operating at energies y/s > 500 GeV would be an ideal place 
to perform precision measurements of a Higgs boson [§,0. 

Supersymmetric (SUSY) theories, in particular the Minimal Supersymmetric Stan- 
dard Model (MSSM) |jlO|, are currently considered as the most theoretically well mo- 
tivated extensions of the SM. Such models predict a rich Higgs phenomenology. The 
MSSM contains 2 neutral CP-even scalars {h^ and H^), a CP-odd neutral scalar A^, 
and a pair of charged scalars {H^,H^). High-energy e^e~ colliders offer a very clean 
environment in which to search for Higgs bosons, and the simplest way to produce a CP- 
even scalar is in the Higgsstrahlung process e"'"e~ —^Z*-^ H^Z and the W boson fusion 
process e+e^ — > H^UeVf. |TT|JT^. The CP-odd A^ possesses no tree-level coupling A^VV 



(where V = Z, W^) and the other tree-level diagrams contributing to e^e ^ A^Z are 



proportional to the electron mass and consequently negligible^. Therefore the dominant 
contribution to e+e^ —>■ A^Z is from higher order diagrams. In previous works we calcu- 



lated the 1-loop induced rate in the context of the CP-conserving THDM [jn| and the 



CP-conserving MSSM ||T5[, finding maximum values of order 0.01 fb and 0.1 fb respec- 
tively if tan/? > 2. The Higgsstrahlung, vector boson fusion, and Higgs pair production 
mechanisms (e^e~ -^ A^h^/H^ and e+e^ -^ h^H^), have all been extensively studied 
(including radiative corrections |]T6|-|T8[1) at NLC energies in the context of the MSSM. 
In recent years the phenomenology of the MSSM with complex SUSY parameters has 



received growing attention ||T9[, [|3], |^, |2^. Such phases may allow baryogenesis [^ , 



and do not necessarily violate the stringent bound from the non-observation of Electric 
Dipole Moments (EDMs) ||2^, I^Sf- The presence of SUSY phases induces mixing between 



the CP-even and CP-odd scalars, resulting in the mass eigenstates H^, iJg ^^^ ^s which 
are mixed states of CP. This mixing affects their phenomenology at present and future 
colliders, both in production mechanisms and decay partial widths ||2^. In this paper we 
will study the processes e"'"e^ -^ Z* ^ H^Z and e~^e~ -^ H^uJ^^ in the context of a NLC. 
Both mechanisms are mediated by the tree-level couplings H^VV, but their cross-sections 
have different phase space and ^/s dependence. In the CP conserving case one of these 
couplings would be zero, corresponding to the absence of the coupling A'^VV. We will also 
study the production of neutral Higgs pairs, e^e" —>■ H^H^ {i ^ j). In the CP conserving 
MSSM, only the vertices Zh^A^ and ZH^A^ exist at tree-level while in the CP violating 
scenario all three couplings ZH1H2, ZH'^H'^ and ZiJ^-^s ^ire generated at tree-level. 
Therefore an observable signal for all i{= 1, 2, 3) in a given mechanism {e^e~ -^ ZH^, 
e^e" — > HfHj or e^e~ -^ H^VePe)-, would be a way of probing CP violation in the Higgs 



sector. Such an approach was used in the context of the THDM in |2^, with related 



analyses in p8[. We will calculate the tree-level rates of the above mechanisms in the 



^Note that at a muon collider |13|, the tree-level diagrams for e^e -^ A^Z cannot be discarded 
anymore and higher order diagrams would induce corrections to the tree-level rate. 



context of the CP violating MSSM, showing that in the most favourable scenarios this 
way of probing scalar-pseudoscalar mixing can be effective if M^^.^ < 300 GeV. 

Our work is organized as follows. In section 2 we outline our approach for evaluating 
the cross-sections for the above production mechanisms in the CP violating MSSM. In 
section 3 we present our numerical results and section 4 contains our conclusions. 

2. SCALAR-PSEUDOSCALAR MIXING IN THE MSSM 

The tree-level Higgs potential of the MSSM conserves CP, which ensures that the 
three neutral Higgs eigenstates can be divided into the CP-even h^ and H^ and CP-odd 
A^. Recent studies ||l9| , |20|p^ have shown that the 1-loop effective potential may violate 



CP resulting in three Higgs mass eigenstates which cannot be assigned a definite CP 
quantum number, denoted by if °, if 2 ^^^ H^ (in ascending order of mass). In the above 
studies it is shown that the CP violation is generated by complex phases which reside in 
the /i term and the soft SUSY breaking parameters At and A^y. These phases generate 
terms M.\p in the 3x3 neutral Higgs mass squared matrix M.f, which mix the CP-odd 
and CP-even scalar fields. These may be given approximately by ||19 

lAP l/ip sin20c'p|A||/i| 






6, 



^lc/5y'tan/?M|j;cjy' sin^cpAfJ^^y 



(4.1) 



where (pcp = arg(y4(/i), and we have only displayed the contributions from the top squarks, 
ti_2, which are dominant for small tan/5. Sizeable scalar-pseudoscalar mixing is possible 
for large |/i|,|74(| > Msusy- In Pll the mass matrix Jvif- is evaluated to one-loop or- 
der using the effective potential techniques and includes large two-loop non-logarithmic 
corrections induced by one-loop threshold effects on the top and bottom quark Yukawa 
coupling. The public code which we will employ in our numerical analysis can be found 
in [^]. If the SUSY phases are set to zero, the mass eigenstates become definite eigen- 
states of CP. A phenomenological consequence of the scalar-pseudoscalar mixing is that 
all the eigenstates H^ possess a tree-level coupling VVH^ with respective strength Cj 
(normalized to SM strength). These couplings can be easily obtained from the covariant 



derivative of the Higgs fields in which the non-physical Higgs fields are expressed in terms 



of the physical mass eigenstates. Following the convention of |22|, the Cj are given by: 

C, = Oi,cos/5 + 02^sin/5 (4.2) 

Here Oij is the orthogonal matrix which diagonalizes M-lj- One can easily show that 

Cl + Cl + Cl = l (4.3) 

Note that this sum rule applies to the tree-level vertices, although it holds to a very good 
approximation if higher order corrections to the vertices are included. We will present 
results using the tree-level values for Ci which we will generate by use of the program 
cph.f [H . 



In the CP conserving MSSM scalar-pseudoscalar mixing is absent. In this case Ci = 
sin(/3 — a), while one of 6*2,6*3 is identified as cos(/? — a), and the other is identically zero. 
C*j = corresponds to the couplings A^ZZ and A^WW, which will take on a non-zero 
value at 1-loop. Hence to lowest order in the MSSM only Hf and one of H2, H^ can be 
produced in the Higgsstrahlung and WW fusion mechanisms, e^e" ^ Z* ^ HfZ and 
e+e~ -^ H^VeT^e- In the presence of SUSY phases all H^ may be produced at tree-level via 
e+e~ — > Z* — > H^Z, and an observable signal for all three H^ would be evidence for CP 
violation in the Higgs sector. We stress here that the smallest of o"(e+e^ ^ Z* ^ ^i^) 
should exceed the maximum rate for e~^e^ — > A^Z in the context of the CP conserving 
MSSM |T^, since the latter would constitute a "background" to any interpretation of 
scalar-pseudoscalar mixing. Note that the process e~^e~ -^ H^uJ^e proceeds via the same 
couplings C*j, but possesses a different phase space and ^/s dependence. This mechanism 
is competitive with the Higgsstrahlung process, and becomes the dominant one as y/s 
increases. We will also consider the mechanism e"'"e" -^ HjH^, which proceeds via the 
coupling Cjk, where C*j = Cjk for i ^ j ^ k. 

In the MSSM (with or without SUSY phases), the properties of the lightest eigenstate 
H^ become very similar to that of the SM Higgs boson in the decoupling region of Mh± > 
200 GeV. In this region, C*i is very close to 1, and so the sum C2 + C^ is constrained 
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to be small (Eq. ([4.3|) ). Therefore we expect that Mh± < 200 GeV will allow larger 
values for the sum C2 + C^, and thus observable rates for both if 2 and if 3 in the above 
mechanisms. Following the approach of ||28| we will take the threshold of observability as 
Cobs = 0.1 fb. This would give 50 raw events (before cuts) for the assumed luminosities of 
500 fh^^. Distinct signals for all three H^ in a given channel would be evidence for scalar- 
pseudoscalar mixing. Note that ]15|] found maximum values of o"(e~*"e~ -^ A^Z) = 0.1 fb 



in the context of the CP conserving MSSM. 

A caveat here is that extended Higgs sectors with more than two doublets or extra 
Higgs singlets (e.g. the NMSSM) would also predict multiple signals in these mechanisms 



2^. We will show that the CP violating MSSM can only produce multiple signals below 
a certain mass for Mh2 and Mh^, and thus any such signal for a larger Higgs mass would 
be evidence against the CP violating MSSM. Therefore the measurement of the Higgs 
mass may act as a discriminator among the models. 

3. NUMERICAL RESULTS AND DISCUSSION 



We now present our numerical results which we will generate with the fortran program 
cph.f [^. We note that this program does not include x^^^^H^ contributions to the 
1-loop neutral Higgs mass matrix, which have been shown to be sizeable in some regions 
of parameter space 0] . The analytic expressions for the various cross-sections are given 
in the literature [Q . For the WW fusion process we will use the exact expression given in 

Graphs showing the numerical values of the couplings Ci and Cij have appeared in Refs 



[|^,Q . However, these papers were more concerned with light Mh2 and Mh^ of interest at 
LEP2 and the Tevatron. Detection at these colliders would require quite sizeable values 
for C| and C|. We are concerned with a NLC collider which has the ability to probe 
cr(e+e~ -^ HfZjHfHjjHfueU) > 0.1 fb, and so we are also interested in smaller values 
for C^iC^ and larger Mh^- An earlier discussion of the potential of a NLC to probe 



scalar-pseudoscalar mixing can be found in [^. However this study only addressed the 



pair production process e^e~ -^ H^Hj and numerical results were only presented for the 
couplings ZHfH^. We shall be presenting results for the production cross-sections of 
all the above mechanisms, in contrast to p9| , |2T|j2^ which were more concerned with the 
numerical values of Ci and Cij and LEP2 phenomenology. Since the WW fusion process 
is considerably more important at NLC energies than at LEP2 energies, our analysis is 
complementary to that in |^l| and extends that of |^2|. For alternative ways of probing 



CP violating SUSY phases at e+e colliders, see p3 |. 



As noted in the introduction, the cross-sections for the processes e^e -^ Zh^/H^, 
e+e- -^ Ah^/H^, e+e" -^ /i°if°, e+e" -^ ZA^ and e+e" -^ VeVeh^/H^ in the CP 



conserving MSSM are accurately known [|14| , p!6| -[lql . Deviations from these rates would be 
evidence for scalar-pseudoscalar mixing. 

The presence of large SUSY phases can give contributions to the EDM which exceed 
the experimental upper bound. To avoid conflict with experiment one may assume that 
the masses of first two generation of squarks are well above the TeV scale while the third 
generation may be relatively light (< 1 TeV) ||2^. A recent paper |24[ suggests that 



sizeable scalar-pseudoscalar mixing would prefer the cancellation mechanism [^ over 
the above mechanism as the solution to keep the value of EDM within the experimental 
limits. Another option is to adopt a non-universal scenario for the tri-linear couplings Aj 



35|] . In particular, one may require arg(/i) < 10~ and Af = (0,0, 1)^4, with At,Ai, and 
Ar taking maximal phases. Such a scenario comfortably satisfies the EDM constraints. 
Since the scalar-pseudoscalar mixing ~ (pcp = argl^t/^), it is sufficient to have maximal 



phase in At to maximize (pcp- However, two-loop Barr-Zee type diagrams [Q can violate 



the EDM constraints for large tan/3 (> 30). Therefore we will restrict ourselves to low to 
intermediate values of tan/3. 

In our numerical analysis we will choose the CP violating benchmark scenario (CPX) 
which was introduced in [^ and maximizes the CP violating effects. The CPX scenario 



is as follows: 

MQ = Mt = Mb = MsusY = 0.5 ^ ITeV , /i = AMsusy 



\At\ = \Ab\ = 2MsusY , \m~\ = ITel/ and |mg| = \m^\ = 0.3Te1/ (3.1) 

Note that fi will be taken real while we allow a CP phase in the soft tri-linear parameters 
At and Af, and in m-. The CP phases of At and Af, are chosen to be equal and may be 
maximal. In addition, we choose the charged Higgs mass and tan/? as free parameters. 

Our strategy to probe the scalar-pseudoscalar mixing requires the identification of 
the Higgs signals as distinct resonances. The inclusion of the phases in At and Ab breaks 



the near degeneracy among Mh2 and Mh^ ||22|, and gives sufficient splittings to allow 



identification of separate resonances for H2 and H^. These splittings may be > 10 GeV, 
which is sufficiently large for a NLC f^ to resolve the separate peaks. This will lead to 
three different peaks in the Higgsstrahlung and WW fusion processes and motivates us to 
present the individual cross-sections for e'^e" -^ ZH^ and e~^e~ -^ vVeH^ for 2 = 1,2, 3. 



This is in contrast to [|3^ where the study was devoted to LEPII energies and the cross- 



sections were summed over the three Higgs states. It has been shown in p!9|j2T| , |22[] that 
the inclusion of SUSY phases may drastically change the size of the couplings ZZH^ and 
ZH^H2 for low and intermediate tan/3. In such cases the bound on the light Higgs boson 
obtained at LEPII may be weakened to < 60 GeV for large CP violation in the MSSM 
Higgs sector. We study the potential of a NLC to discover such a weakly coupled Higgs. 
In Fig. 1 the left (right) plots depict regions of a{e'^e~ — > ZH^) in the plane 
{MH,,a.ig{At)) for ^s = 500 GeV, tan/3 = 6(15), and Msusy = 1000(500) GeV. In 
all plots the charged Higgs mass has been varied in increments from 140 — > 400 GeV, 
which determines the values of Mh^- Comparing the left and right plots it is clear that 
lower tan /3 provides larger cross-sections for H2 and H^ over a wider region of the plane, 
corresponding to the fact that the scalar-pseudoscalar mixing is enhanced. For H^ dis- 
covery is possible over most of the {MHi,Si.Tg{At)) plane, with small unobservable regions 
where a{e^e^ -^ Z^i) < 0.1 fb which occur for arg(y4() ^ 1-5(2) for tan/3 = 6(15) and 
Mhi ^ 105(115) GeV. The smallness of cr(e+e^ -^ ZH^) is due to the suppression of Ci. 
This suppression arises when O21 changes sign, which induces destructive interference in 
Eq. (|4^ ). Note that Ci is dominated by O21 and so Ci also flips sign for tan/3 = 6(15) 



and arg(Ai) ^ 1.5 (~ 2). C2 is positive in both cases tan/? = 6 and tan/? = 15, and is 
maximized for M^^ < 150 GeV; it is minimized for aig{At) ~ l-5(> 2) and Mh2 ^ 150 
GeV for tan/? = 6(15). For H2 and iJ^, both o"(e^e~ -^ Zi/23) can be observable over 
a wide region of the plane, even up to relatively large mass values e.g. for tan/? = 6 and 
arg(A) = 1, (r{e+e- -^ ZH^^) > 0.1 fb for Mh., < 250 GeV and Mh,, < 270 GeV. Note 
that the scalar-pseudoscalar composition of H2 and H^ can change with increasing Mh- 
e.g. for the tan/? = 6 plot with low arg(A() (i.e. small scalar-pseudoscalar mixing) one 
can see that H^ is dominantly scalar for low masses, and has a much larger cross-section 
than for that for H2- As Mh2 3 increases, H2 has the larger scalar component and may be 
produced with an observable rate for M^^ < 300 GeV. The coverage for e^e~ — > HfueV^ 
at the same y/s is comparable to that in Fig.l. 

Fig. 2 shows a{e+e- -^ ZHf) for v^ = 800 GeV, tan/? = 6, and Msusy = 1000 GeV. 
Ones sees that the coverage is inferior to that in Fig.l since a{e^e~ -^ ZHf) is reduced 
for larger \/s. 

Fig. 3 is analogous to Fig.l and shows (j{e^e^ -^ Hfve^e) for y^ = 800 GeV. Here 
we find improved coverage compared to that in Figs. 1 and 2, since the cross-section for 
this process is enhanced with increasing ^/s. For the tan/? = 6 plot with arg(Af) = 0.5, 
both a{e^e^ -^ i/^s'^e^e) > 0.1 fb for M//23 < 300 GeV. The window of unobservability 
for H^ has essentially been closed. 

In Fig. 4 we show a{e^e^ -^ H^H^) in the plane (M^., arg(Ai)) for tan/? = 6, 
Msusy = 1000 GeV, and ^/s = 500 GeV. This mechanism offers comparable cross- 
sections to those for the Higgsstrahlung and WW fusion processes, and consequently is 
also effective at probing scalar-pseudoscalar mixing. The best coverage is obtained at 
arg(At) ^ 0.5. Using the fact that Cf = C'^f, for i ^ j ^ k, the behaviour of pair 
production e"'"e' -^ H^Hj can be roughly understood from the rate of the Higgsstrahlung 
process e+e^ -^ ZH^. As can be seen from the plots, there are some similarities between 
e+e" -^ H^H^ and e+e" -^ ZH^, for i ^ j ^ k. Note that a{e+e- -^ H^H^^) have large 
cross-sections (> 5 fb) in the region where the Higgsstrahlung and WW fusion processes 
have very suppressed rates. This situation corresponds to strong scalar-pseudoscalar 
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mixing, since in the absence of SUSY phases one of a(e+e -^ H^H2^ would be zero. 

4. CONCLUSIONS 

We have studied the production processes e'^e^ -^ H^Z, H^Hj and H^UeV^, in the 
context of the CP violating MSSM. We showed that in a given channel the cross-section 
for all H^ {i = 1,2,3) can be observable at a Next Linear Collider and would provide 
evidence for scalar-pseudoscalar mixing. At y/s = 500 GeV the coverage of e+e^ -^ H^Z 
and H^UgVe are comparable, with observable cross-sections for Mh2 ^ 250 GeV and 
M//3 < 270 GeV for the most favourable choice of aig{At). At ^/s = 800 GeV, the process 
e~^e~ -^ H^UeVe offers superior coverage, with a reach up to Mh2 3 < 300 GeV in the most 
favourable cases. The scalar-pseudoscalar mixing causes a mass splitting between H2 and 
H^ which should be sufficient for separate peaks to be resolved at a NLC. The problematic 
region of a light H^ with a very suppressed coupling to vector bosons (VVHi) has a 
window of unobservability at ^/s = 500 GeV and arg(A() ^ 71 /2. This is almost closed 
at ^/s = 800 GeV in the HfueUe channel. The mechanism e^e~ -^ H^H^ is competitive 
with the above mechanisms for probing scalar-pseudoscalar mixing at ^/s = 500 GeV, 
and can comfortably detect Hf in the region of suppressed coupling VVH^. 
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